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Constant- Weight and Constant-Charge Binary 
Run-Length Limited Codes 

Oleg Kurmaev, Member, IEEE, 
Abstract 

Constant-weight and constant-charge binary sequences with constrained run length of zeros are introduced. For these sequences, 
the weight and the charge distribution are found. Then, recurrent and direct formulas for calculating the number of these sequences 
are obtained. With considering these numbers of constant-weight and constant-charge RLL sequences as coefficients of convergent 
power series, generating functions are derived. The fact, that generating function for enumerating constant-charge RLL sequences 
does not have a closed form, is proved. Implementation of encoding and decoding procedures using Cover's enumerative scheme 
is shown. On the base of obtained results, some examples, such as enumeration of running-digital-sum (RDS) constrained RLL 
sequences or peak-shifts control capability are also provided. 

L Introduction 

The sequences with constrained run length of zeros are known in literature as dk sequences. In these sequences, single ones 
are separated by at least d, but not more then k zeros. A dkr sequence is a dk sequence, ending in a run of not more then r 
trailing zeros. A dklr sequence is a dkr sequence, beginning with a run of not more then / leading zeros. These sequences, 
their properties, and applications are described in [1] in detail or, briefly, in comprehensive overview paper [2]. 

A general enumerative scheme for encoding and decoding binary sequences has been presented by Cover [3]. We use this 
technique for determining the number of constrained sequences. 

Let {0, 1}" be the set of all binary sequences of length n and let x — {xi,X2, ■ ■ ■ , Xn) denote a generic element of this set. 
Let S''{n) = {a; e {0, 1}" | satisfies the d, k, r constraints} and let ^{n) = {a; G {0, 1}" | satisfies the d, fc, I, r constraints}. 

Using Cover's method, cardinality of ,S^{n) can be computed as shown in [4]. To do this, the number of dkr sequences, 
which begin with one, is calculated as 

fc+i 

^ \,9'{n-j)l n>d+k. 
j=d+i 

Then the number of dklr sequences is calculated as 

min(n,^) 

By = J2j=i denote the weight of the sequence x. The number of unconstrained constant-weight sequences may be 
simply obtained as ("), see [5]. Methods for calculating the number of constant-weight dkr sequences is given by [6], [7]. 
Under NRZI encoding [2] we understand mapping the source sequence x to bipolar sequence z, z E { — 1, 1}" such that 

_ Jzj-i, = 0, 

^ ~ I _ — 1 

[ Zj-ii — 1, 

zo = 1. 

By cr = X]j=i Zj denote the digital sum or charge of the sequence z. Observe that a G [—n,n], where cr admits even 
values whenever n is even and odd values whenever n is odd. For example. Table |T] shows sequences x and z, theirs v and a 
respectively. 

Our central goal will be to determine the number of constant-weight and constant-charge run length limited sequences as 
well as an exact estimation of this number We also intend to show that using Cover's enumerative method provides us with 
necessary values for error control. Namely, these values follow from weight distribution, run-length distribution, and charge 
distribution, which may be obtained by Cover's technique. 

Although we will consider constant-weight and constant-charge RLL sequences together, our main efforts will be focused on 
constant-charge RLL sequences. Some results concerning the constant-weight RLL sequences are known with the contributions 

The material in Section|ll]of this paper was presented in part at the 10th International Workshop on Algebraic and Combinatorial Coding Theory (ACCT-10), 
Zvenigorod, Russia, September 2006. 

The author is with the Moscow Institute of Electronic Engineering (MIEE/MIET), 124498, Moscow, Russia, (e-mail: <kurmaev@org . miet . ru>) 
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TABLE I 

All Lexicographically Ordered dklr 



Sequences of Length n = 8, With 





Constraints d = 


2, k 


= 4, 


/ = 


1, r = 


3. 






X\ 




^3 




X5 




X7 


xs 


u 






21 


22 


23 


24 


25 


26 


27 


28 




cr 








1 














1 





2 






1 


-1 


-1 


-1 


-1 


-1 


1 


1 




-2 




1 





1 











1 










2 







1 


-1 


-1 


-1 


-1 


1 


1 


1 




2 





1 








1 











2 






1 


-1 


-1 


-1 


1 


X 




I 




2 


3 





1 








1 








1 


3 







1 


-1 


-1 


-1 


1 


1 


1 


-1 




4 


1 














1 








2 






-1 


-1 


-1 


-1 


-1 


1 


1 


1 




-2 


5 


1 











1 











2 







-1 


-1 


-1 


-1 


1 


1 


1 


1 




6 


1 











1 








1 


3 






-1 


-1 


-1 


-1 


1 


1 


1 


-1 




-2 


7 


1 








1 











1 


3 







-1 


-1 


-1 


1 


1 


1 


1 


-1 




8 


1 








1 








1 





3 






-1 


-1 


-1 


1 


1 


1 


-1 


-1 




-2 



" By Af we denote the lexicographic index of se- 



quence. 



coming from Lee [8], Forsberg and Blake [9], Ytrehus [6]. We will cite some of their results for the sake of generahzation 
and in comparison with the similar results for constant-charge codes. 

Run-length distribution does not considered in this paper, although the length of each run is accounted in Cover's enumerative 
technique. The problems of defining this distribution and accompanying problems are studied in source coding. The reader can 
find these materials in various publications since weU-known Huffman's paper [10]. 

The rest of this paper will be organized as follows. First, in Section [III we derive recursion relations for determining the 
number of the sequences. Further, in Section |llll we obtain direct formulas for the same. Next, in Section IIVI we consider 
the numbers of our sequences as coefficients of formal power series and derive generating functions. We also prove that 
generating function for enumerating constant-charge RLL sequences does not have a closed form. Then, in Section |V] we give 
an enumerative algorithm for encoding and decoding these sequences. Some remarks and application notes are presented in 
Section |VT] Finally, conclusions are drawn in Section I VII I 

II. The Number of Sequences 

Consider run-length constrained binary sequences of length n and weight v. Let A^^ be the number of these sequences which 
begin with one. Let Aj^ be the number of these sequences which begin with a leading run of zeros. 

Since an internal run of zeros succeeds a leading run of zeros then the leading constraint I does not affect A^. For convenience, 
below under A'^ we consider A'^{d, k, r) and under AJ^ we similarly consider A'^{d, k, I, r). 

Suppose that a unique sequence of zero length and zero weight exists. Let it be a sequence which begins with one. Then 

Al = l. 

Proposition 1. The numbers A'^ and A'^ can be obtained as: 
If V = 1 and the sequences begin with one, then 

10, otherwise. 

If V > 1 and the sequences begin with one, then 



If V = and a leading series of zeros is running, then 



^JO, n<d+l, 



0, otherwise. 



If v > Q and a leading series of zeros is running, then 



min(n,^) 



3 



TABLE II 
An Example of Charge Changing 



(a) After a leading one. 



(b) After leading zeros. 



m 


X2 Xs, X4 X5 


0"m 


5 


10 10 
-1 -1-11 1 


-1 


n 


^1 X2 X4 Xg X7 Xg 


0"n 


8 


1 1 1 
-1-1-11 1 1 -1 -1 


-2 



m 








Xl 


X2 


2:3 


X4 


^5 


0"m 


5 








1 








1 















-1 


-1 


-1 


1 


1 


-1 


n 


Xl 


X2 




X4 


^5 


XQ 


X7 


a;8 




8 











1 








1 









1 


1 


1 


-1 


-1 


-1 


1 


1 


2 



Here d > 0, k > d, I > 0, r > 0. Proof: In the case of v = 1, there is only a traihng run of zeros in the sequence. It 
gives us the only allowed sequence which length lies in the interval [1, r + 1]. Therefore, ([T]l is evident. In the case of > 1, 
according to Cover's enumerative method [3], we build the recursion by the following way. Let us consider a possible run 
of zeros, which follows the leading one, as a prefix for the following subsequences beginning also with one. Assuming the 
length of the prefix grows from d + 1 to min(n, fc + 1) and weight of this prefix equals one, we obtain the number of these 
subsequences evidently equal A'^^z] for each allowed prefix and the total number is expressed by 

The other case is when a leading series of zeros is running. In the case of zero weight, the leading run of zeros is the trailing 
one. This also gives us the only allowed sequence which length lies in the interval [0, min(/, r)]. Hence, (O is also evident. In 
the case of nonzero weight, there exist only zero weight prefixes which length lies in the interval [0, min(rt, I)]. Subsequences 
beginning with one follow this prefixes, thus, the total number of sequences is expressed by (01). 

Consider bipolar run-length constrained sequences of length n and charge cr. We can do this in terms of source sequence 
X. This allows us to obtain results similar to Proposition [T] In this case a = where i^j ~ ^i- Let be 

the number of these sequences, which begin with one. Let C!^^ be the number of these sequences, which begin with a leading 
run of zeros. 

Since an internal run of zeros succeeds a leading run of zeros then the leading constraint I does not affect C,j . For convenience, 
below under C,'^ we consider C!^^(d, k, r) and under C*^ we similarly consider C!^{d, k,l,r). 

Proposition 2. The numbers C"-^ and can be obtained as: 
If a = —n and the sequences begin with one, then 

= + (5) 

I 0, otherwise. 

If a ^ —n and the sequences begin with one, then 

^/O, n<d+l, 



If a — n and a leading series is running, then 



If <T ^ n and a leading series is running, then 



1, n < min(^, r), 
0, otherwise. 



inin(n,/) 

C= E ^rj- (7) 

i=o 

Here d > 0, fc > d, Z > 0, r > 0. Proof: Is similar to Proposition [T] except the differences between the weight and 
charge. Namely, if the sequences of length n begin with one, then the charge an can be obtained as 

where m and am is the length and charge of following subsequence that also begins with one. Since m = n ~ j, therefore, 
the number of subsequences in (|6]l must be C^'^J-' . 

If the sequences of length n begin with zero, then the charge (T„ can be obtained as 

(T„ = n - m + (T,„. 

Thus, the number of subsequences in (|7]l must be C^Z] ■ For example, the charge changing is shown in Table |II] 
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TABLE III 

An Example of Weight and Charge Distribution 



(a) For sequences beginning with one. 
Constraints: d = 2, fc = 4, r = 3. 



(b) For sequences beginning with a leading run of zeros. 
Constraints: d = 2, k = i, I = 1, r = 3. 







'-'n 


u, a 


012345678 


-8 -7 -6 
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-3 


-2 


-1 
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7 8 


n 
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1 
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1 
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10 








1 
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110 






1 




1 
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2 









1 




1 
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3 
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1 


1 
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3 10 










1 
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1 
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0023 00000 












3 


2 



























v, a 


012345678 


-8 -7 -6 


-5 


-4 


-3 


-2 


-1 


1 
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6 7 8 


n 































1 












1 
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1 1 
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2 
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1 
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2 








1 




1 
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2 10 






1 




2 
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13 









2 




2 
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5 
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2 


2 
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6 10 










2 
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2 
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0054000 












4 


4 




1 













Using relations derived in this section, we can write the weight and charge distribution of our sequences. It seems to us that 
mostly convenient form for presenting this distribution is a triangle table, like Pascal's triangle. An example of such distribution 
is shown in Table Hill 

We may consider (|6]l as an implicit mutual recursion. Indeed, an element in left slanting row in the triangle of charge 
distribution (see Table Hllf a)) depends on elements in right slanting rows and vice versa. In more details, this will be shown 
in the next section. 



III. Direct Equations for the Number of Sequences 

Consider dkr - limited sequences. The problem of defining the number of constant-weight sequences, in which term's 
coefficients are Fibonacci numbers, i.e., for d - sequences, has been solved by Riordan [11]. He presented a direct, not recursion 
method for calculating the number of constant-weight d - sequences. Lee in [8] and Ytrehus in [6] extended Riordan's method 
for the dkr - constrained sequences. Our method of deriving the similar direct equations rests on recursions from the previous 
section. 



A. Calculating the number of constant-weight sequences 
From ([T]i we have that 

(suppose that (^) = if /c > n). From recursion relation (O we have 

— + ^n-(d+2) + ' ' ' + A\-(fe+l)' 

= ^ri-(ci+l) + ^n-(d+2) + ^ ^n-(fe+l)> 



— ^n-(d+l) ^ n-(d+2) ^ ^ n-(fe+l)- 



Substituting A\ into A^, we obtain 



^2 _ ^ri-(rf+l)-lj _ ^n-(d+l)-l-(r+l)'j 

_|_ ^n-(<i+2)-lj _ ^n-((i+2)-l-(r+l)-j 

^ ^n-(fc+l)-lj _ ^n-(fc+l)-l-(r+l)-j 

_ (n-{d+l)-\^ _ ^n-(r+l)-(ci+l)^ 

_ ^|^n-l-(/c+l)j _ ^ri-l-(r+l)-(fe+l)-j^ 



5 



Likewise, 



^3 _ ^n+l-2(d+l)j _ ^n+l-(r+l)-2(d+l)^ 



n-(d+l)-(fc+l)-j |-ji-(r + l)-(d+l)-(fe+l)N 



?i-l-2(fe+l)^j _ ^ri-l-(r+l)-2(fe+l)-j 



^j; = E(-i)^(Y)(r 



i+(i.-i-j)-(i.-i-j)(£;+i)-i(fe+i)\ 



_ /n-l + (l/-l-i)-(r+l)-(l/-l-J)(d+l)-i(fc+l)^ 

1/ > 1. (8) 

By q denote the number of possible lengths of runs beginning with one 

q = k- d+l. 

Then we can rewrite ^ in more convenient form 



J=0 



(n-l-(r+lM.-l)<i-J9)\ J.>1. (9) 



B. Calculating the number of constant-charge sequences 

Now, consider the constant-charge sequences. For these sequences we may simply show that 

This directly follows from initial conditions (|5]l of recursion relation 

Thus, we define the first left slanting row in our triangle of numbers of constant-charge sequences. For example, see 
Table HHa). 

We also may show that for right slanting rows 

C° = 1, Cl =0, = 0, . . . d > 0, 

': : : : (11) 

c-;r = 1, c;^2^'^^' ^ 0, c;l;^'^^' = o, . . . d > r + 1, 

^-r+2 _ r, ^-(r+2) + l _ „ ^-(r+2)+2 _ „ 

Let S be an even number which shows how the charge a differs from n or from —n. In other words, S/2 is an index of a 
slanting row in the triangle of numbers of constant-charge sequences. 
Therefore, we can rewrite ( fTOl l as 

^.^+s_ffS/2\ fS/2-1 



V / V 



n - 1 - (r + 1) 



_oy V 

Also we can rewrite (fTTT i as 

^„_,_ ^^<5/2-l^ ^(5/2-l-(r+l) 



(12) 



V y V 

j ~ I 
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If d < (5/2 then for right slanting rows, we have from recursion relation (|6]l that 

f^n-S _ f^-n+S-(d+l) ^-n+S-(d+2) ^-n+5-(fc+l) 
— + '^n-{d+2) ^ ^ ^ri-(fe+l) 

_ ^-(n-(d+l))+S-2(d+l) ^-(n-(d+2))+(5-2(d+2) 

~ ^ri-(d+l) + *^n-(d+2) 

This confirms a mutual nature of (|6]l. The series ( fT4] i terminates early if cr < — n, i.e., whenever someone of (5 — 2((i + 1), 
5 — 2{d + 2), . . . , S — 2{k + 1) becomes less then 0. Therefore, combining (fTsT i with (fl4] i. we obtain 

rin-S _ ^-(«-(d+l))+<5-2(d+l) ^-(ri-(d+2))+(5-2(d+2) 
~ +^ri-(d+2) 
^-(n-(fe+l))+5-2(fe+l) 

+ ((^^r^) - c^'-V'^'^)) (r-/^) - ^)) 

_ ^|^(5/2-(d+l)j _ |-(S/2-l-(<i+l)^ 
_|_ ^A72-(d+2)j _ ^<S/2-l-(£i+2)j 
|_ ^''/a-^fc+l)^ _ ^5/2-l-(fc+l)j^ 

X (("-^/^) - r-V-^)) , <5>0, 



then 



™_a _ ^-(n-(d+l))+5-2(d+l) ^-(n-(d+2))+5-2(d+2) 
'-'n — '-'n-(d+l) +^n-(d+2) 



(( 



^-(n-(fc+l))+<5-2(fe+l) 



X (ro^') - r'o^'"')) . '5>0. (15) 

Consider an example. For 6 = 2 and for a = n — 6, we have from recursion relation ^ that 



or 



^n-2 _ ^-ri+2-(d+l) ^-n+2-(d+2) ^-ri+2-(fe+l) 
^ ^n-{d+l) +"-'ri-((i+2) + ' " + "-"n- (fc+1) 

_ ^-(n-(d+l))+2-2(d+l) ^-(n-(d+2))+2-2(d+2) 
~ + ^n-{d+2) 

j_ x^-(n-(fe+l))+2-2(fc+l) 
'-"ri-(fe+l) 



Since this equation specifies the right slanting rows, then we have to substitute ( fTOl i for C„i:" if c? = 0, and ( fTST l otherwise. 
Therefore we obtain the next result 




d = 0, 
otherwise. 



Using (flST l. we can rewrite this piecewise-continuous equation as 



( 



^l-(fe+l)j _ J'l-l-(A:+l)^ 



1 



For 6 = 2 and for tj = — n + (5, we have from the same recursion relation 

^-n+2 _ wn-(d+l)-2 ^n-(d+2)-2 ^ri-(fc+l)-2 
^ "-'n-(d+l) +^n-(d+2) '-'n- (fc+l) 
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and obtain 



X (^fn-l + l-(d+l)^ _ ^i-l-(r+l)-(d+l)^ 



n-l-(fe+l)'j _ ^n-l-l-(r+l)-(fc+l)^ 



^n-l-(fc+l)-j _ ^ri-l-l-(fe+l)^ 

Let A be a value which has the similar combinatorial meaning as A in (|2]) and let 



^r = E(-i)'(7) (( 



^ ^ (S/2+va-j-(m-j)(d+l)-j{k+l)'^ 
^5/2+m-j-l-(m-j)(d+l)-j(fc+l)-j\ ^jg-j 



and 



^r = E(-i)'(7)((' 



^ ^ ^5/2+m-i-l-(m-i)(d+l)-i(/c+l)^, 

_ M/2+m-j-(d+l)-(m-j)(d+l)-j(fc+l)\ 

V m / 

, M/2+m-J-l-(fc+l)-(m-j)(d+l)-j(fe+l)\ 

V m / 

_ ^5/2+m-j-l-(r+l)-(m-j)(<i+l)-j(fe+l)-j^ ^ ^jy^ 



Also, let 



j=0 

_ ^n-5/2+m-j-l-(r+l)-(m-j)(d+l)-i(fe+l)-j^ ^-jg^ 

m 

^n>=E(-i)'(7)(r 



and 



-5/2+m-j-(m-j)(<i+l)-j(fc+l)N| 

-^^ V j" / 

J=0 



_ /ri-i5/2+m-j-l-(m-j)(d+l)-j(fc+l)^ 
V m / 

Using this notation, we may rewrite C~" and C,7"^^ from our examples as 

Cji"^" — ^2-^n,2 + ^2^n,2- 

Likewise, 

0„ — /l4i3„ 4 + /i4i3„ 4 + /l4i3„ 4 + /i4i3„ 4. 

Continuing in the same way, we see that 

5/2 



(19) 



= E (^"^r'M + ^r"'^n>) ■ (20) 
m=0 
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Observe that the upper limit of summation might be written as a quotient of S and 2{d+l). In such case 



-n+S _ 



(21) 



Note that summing C„"+ in right slanting rows does not affect B™^ and B™g. This affects only A™ and A™. Indeed, 
from ( fTSl l it follows that index of summation, which subtracts from n and from 6/2, annihilates in ( fTSl l and fT% and remains 
in (O and ([T7]i. 

Now, consider We have 

k-\-l m 

i=d+l j=0 

^ |^|-5/2+m-i-(m-j)(d+l)-j(fe+l)-(d+l)^ 

_ ,s/2+m-]-l-{m-]){d+l)-jik+l)^id+l)\ 
Km) 
, M/2+m-j-(m-i)(d+l)-j(fc+l)-(d+2)\ 

V m / 

_ M/2+m-j-l-(m-j)(£i+l)-i(fe+l)-(d+2)\ 

V m / 
I tS/2+m-j-{-m-j)(d+l)~j(k+l)~(k+l)\ 

\ ra 
lS/2+m-3-l-(m-3)(d+l)-j(k+l)-(k+lY 
\ m / 

m\ ( (S/2+m-j-(m-j)(d+l)~jik+l) + l-id+l)\ 
rn+1 ) 



3=0 



5/2+r,i-]-{r,i-j)(d+l)-]{k+l)-{k+l)\ 
m+1 / 

(S/2+m-i-l-(m-i)(d+l)-j(A;+l) + l-(d+l)\ 
rn+l / 

/S/2+m-i-l-{m-j){d+l)-j(k+l)-(k+lY 
\ m+1 , 



This means that 



Similarly, 



m\ / M/24-m+l-j-(m+l-j)(c?+l)-j(A;+l)>, 



rn+l 



J=0 



/5/2+m+l-j-l-(m+l-j)(^^+l)-J(fc+l)^ 
V m+1 , 



m+1 



E(-i)^(A)(( 



5/2+m+l-j-(m+l-j)(d+l)-j(fc+l)\ 
m+1 / 



/(S/2+m+l-i-l-(m+l-i)(d+l)-j(A;+l)^ 
V m+1 , 



m+1 
J=0 



^A72+m+l-j-(m+l-j)(d+l)-j(fe+l)N 
^ m+1 } 



/5/2+m+l-j-l-(m+l-j)(d+l)-j(fc+l)^ 
V m+1 , 



k+l 

E 

i=d+l 



k+l 



i=d+l 



Using this result and using (flSl l again we see that 



(5/2 



1-5 



m=0 



(22) 
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and like (1211 1. we may rewrite 

m=0 

For (T = — + (5, we have from recursion relation ^ that 

- ^n-{d+l) + ^n-id+2) ^ ^ ^n-(fe+l) ' ("^^^ 

Similarly, we can see that summing C^~^ in left slanting rows does not affect A™ and A™. This affects i?™^- and i?™^-. 
Namely, from ( l23b it follows that index of summation subtracts only from n. 
Consider X^iil+i ^"i-i s- ^PPlying a reasoning chain like ( |22] | yields 

fc+l jn+l 

E klm = E(-i)'(T) 



X 



-<5/2+m+l-j-(m+l-i)(d+l)-j(fe+l)\ 
m+l J 



This means that 



Similarly, 



fn-S/2+m+l~j-l-{r+l)-{m+l~j)(d+l)-j{k+l)-' 
\ m+l ; 



fc+1 



Ertm rjm+1 



fe+1 



Enm nm+l 

i=d+l 



Using this result and using ( |23T l we obtain 



(5/2 

^n"+'= E '^>0. (24) 



m=0 



This proves the inductive hypothesis ( l20l i. 
Like ( |2TI ). we may rewrite ( l24b as 



n — 

m— 

Substituting n + a for (5 in (l2Tl l. we get (l26l l at the bottom of the page. 



2{d+l) 



E EM)^(7) ( 



m— 



E(~l)"'(™) f (^T^^™'*^^'?) — ^^^^-l-('-+l)-md-jq^ 



J=0 



+ 2 



2^-l-(fc+l)-(m+l)d-jq^j _ ^^-l-(r+l)-(m-l)d-jqN, 



(m , _ 



, »— g -I J ■ . 



(26) 
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IV. Generating Functions for Enumerating the Constant- Weight and Constant-Charge Sequences 

From triangle table (see Table Ulll a)) it follows that there exist two types of sequences and C^. The infinite sequences 
of the first type sae A^, A'(, . . . then Cq , C2 , ... or Cf , CJ, .... The finite sequences of the second type are A}^, A\^ . . . , A^ 
and C~", C"""*"^, . . . , C". Kolesnik and Krachkovsky described in [12] a recursive calculation of generating function for 
enumerating constant-weight sequences of the first type. Lee in [8] suggested a direct method for the same. Here we obtain 
the generating functions of the both types. 

A. Generating Function of Sequence Aq, A'^, . . . 

Consider infinite sequence Aq, A^, A2, ■ ■ ■ ■ Define generating function of this sequence as 

00 

n=0 

Substituting (|9]l for A"^ here and changing the order of summation we have that 



00 



fyl-l-il^-l)d-Jq^^^n 



V 



=(ly-l)+l+(l/-l)fl 



^n-l-(u-l)d-{r+l)-jq^_l_n 

= (i^-l) + l+(i/-l)d 
+ (r+l)+jg 



By using the fact that |] {"T)*" = (^''/(l - tT^^) we get 

n— i/+a 



(i-ty 



(27) 



B. Generating Function of Sequence A^, A^, . . . , 
Consider a finite power series 

n 



We can obtain An{y) in closed form by the following way. Consider a formal power series 

OQ 

1^=1 

We can achieve convergence of this series by assuming t arbitrarily small. Rewrite JZTl l as 



A'^it) 



1-t 



1 - t 



Then we have 



1-t 



yt{i ~ r+i) 



r+1 



1-t 



1-t- ytd'+^ + yt^+'^ ■ 



(28) 



By using the fact that dA{t,y)/dt — ^ nAn{y)t'^ ^ we obtain 

n=l 
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From Cauchy's integral representation, we write 

d'^A{t,y) nl 



yr(l-r'-+i) 



27ri 7r (1 - r - yr'^+i + yT'=+2) (r - 



dr. 



Then 



y 



My) = 



1 - r-^+i 



27ri Jp (1 — r — yT'^+'^ + yr'^^^} t'' 
where F is a counterclockwise simple closed contour surrounding the origin of the complex plane, small enough to avoid any 



dr, 



(29) 



(30) 



other poles of integrand in ( |29] l. The application of the residue theorem yields 



fc+2 _ ^ 



m— 1, 



where ti, T2, . . . , Tfc+2 are roots of 1 — t — yr'^'^^ + yr^^'^ and ri = 1. 



C. Generating Function of Sequences Cq , , . . . or Ci , C3 , . . . 

Similarly, consider an infinite sequence Cq , CJ" , , .... In such case, we define generating function as 



)/2 



n=0 
n even 



ri=l 
n odd 



where i G R. 

It is easily shown that (|26] | can be written as a sum of 12 terms with proper signs. 
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Define these terms C(t)^ as 



where 



and 



i,=i 



r ra+<T 
I 2(d+l 



n+fx , 



771 — 

X iE(-i)'(7)('^ 



-p2-md~jq\ 



P2 



0, 7=1,2,7,8, 

1, / = 3,4, 
Pi = — d, / = 5, 6, 

(fc + l)-d, 7 = 9,10, 

(r + l)-d, 7=11,12, 

'0, 7=1,3,5,7,9,11, 
l + (r + l), 7 = 2,4, 
1, 7 = 6,8,10,12, 

'0, 7=1,2,3,4, 
1, 7 = 5,6,7,8,9,10,11,12. 



(31) 



(32) 



(33) 



Without loss of generaUty we will consider a term of (|32] | as and assume that := C(t)^ below in this section. Of 
course, this implies using (l32T i when calculating. 



12 



We can expand the inner product of two series using double sums as follows: 

(—1)^ ('"7'^) -Pl-md-jij 



3=0 



m—n m 



m \ 



= EE(-i)"^"("r)(:) 



u=0 11=0 



X 2 



n-\-a , n—a , 
—pi—ma—uq~fi\ f — — —p2 — md—vq\ 



then interchanging the sums, we get 

oo /m—1.1 m 

where by g{t) we denote the inner sum. 

oo 

9it)- E 



m=0 \ u=0 v—0 



n-\-(T 



—pi —md—uq—i.L\ 



2 — Fi — "'■Lt— f-ty — /-ij ^ — 2~ 



9{t) 



(34) 



(35) 



n—p 
n even 



odd 



where p = a mod 2. 

Further, we will derive the generating function using arguing style of orthogonal polynomials theory [13]. 
Let G(n) = (^-Pi-"^'*-"?-^') (^-P2-m<i-^9-)^(«-p)/2^ m - > 0. The sequence G(n) may content leading zero 
elements. In this case we can not write a ratio between consecutive terms. To prevent this we can use the following rule 



Er/)r/)- 

J=0 



3=A+B 

CXD 

(.J=C+D 

Let a be the A + B — C — D in the conditionals above. In our case a = —a + pi — P2 + (u — v)q and we can rewrite ( |35] ) as 

— /3 + A /m — /i + /3 + |a| + 



oo 

xE 

J=0 



m — [3 



m — p, + 13 



(36) 



where 



and 



[-cr/2j +pi a>0, 
[(t/2J +P2 + 1"?, otherwise 



/9 = 



/i, a > 0, 
0, otherwise. 



Now, consider the sum in 



By G{j) denote a term of this series '^t^i^+p^^^)t^ ™d write the term ratio as 

G[J + ij _ ( m-ff A m-Ai+/3 / 
U + 1) 771 — P 

{\a\ +j + Ijm-p+a 
_ {m- f3 + j + l){m- fi + (3+ \a\ +j + l)t 
~ ij + l){\a\+j + l) ' 
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where {a)j denote the Pochhammer-symbol which is defined by (a)o = 1 and {a)j = a{a + l)(a + 2) ■ ■ ■ {a + j — 1), 
i = 1, 2, 3, . . . Obviously that G(0) = and hence, 

By summing G{j) we get the hypergeometric series as follows: 

g{t) = t'"('^+l)+'^G(0) 



/m-/3+l,m-At + /3+|a| + l ,\ 



If we apply the Pfaff-Euler transformation formula [14], then we obtain the next terminating series 
By several transformations we get 

(m - /? - |a|)(m - /? - |q| - 1) ■ ■ ■ (m - /3 - |a| - j + 1) . 
(|a| + l)(H+2)...(H+l+i-l) 

then, using the distributive law, we obtain 

X (m - M + /? + |a|)(m - + /3 + |a| - 1) • • • (|a| + j + 1) 
X (m - /3 - |a|)(m - /? - |q!| - 1) • • • (m - /3 - |a| - j + 1)*-'. 

Now, our goal is to apply the Leibniz rule to this finite series. We have a power series in one variable t. But, we need a power 
series in two variables because we have rising and falhng factorial powers here. Recall the Mobius transformation and suppose 
that 

X — 1 

t = , (assume x till the end of this section) (37) 

x + 1 



and rewrite 



„U) ^ ^m(d+l)+o-|a| (^_+l2 



m-M+/3+l 



22m-^+l 
m—fj,-\-l3 

E 



m — ^ + P 



{m-fi + py. ^ V j 

x{m- IJL + I3+ \a\){m - 11 + +\a\ - I) ■ ■ ■ 
X (m - /i + /3 + |q!| - (m - + /3 - j) + 1) 

X (a; — i)"*~''+/5+l"l~('"~'"+/'~-?) 
X (m - /3 - |Q;i)(m — 13 - \a\ - 1) ■ ■ ■ 
x{m-p-\a\-j + l){x + 
So, we obtain {m- fi + P- j)th derivative of (a; + ^^^^j^ derivative of (a; + 



X E 



m 



22m-p+l (to-^ + /3)! 

X — ^ , 

dxi 
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TABLE IV 
Nomenclature for |42j 



The first double sum: 


The second double sum: 


Qu = —cr + pi — P2 + wq, 


(41) 


OLv = —a + Pi — P2 — f 9, 


(42) 


a„ = [a 2} +P2 + < 

1 Ou, else, 


(43) 


a. = L-<x/2j+pi-|^^ else. 


(44) 


such that a = wq + a^. 

j /i, a„ > 0, 
'^^ " |0, else, 


(45) 


such that a = wq + a„. 

J a« > 
" |0, else. 


(46) 



where a = a—\o\. Finally, we have the Rodrigues' formula 




where 



(38) 



+ 1/ 4 2(1 -i)' 

Since we suppose ( [37] l. then by r = |^ we denote the similar mapping in the complex plane. From Cauchy's integral 



representation, we write 



1 /■ T-|a|+/3(^2 _ -j^)"!. 

^ ^ (^_l)M(^_a;)m-p+/3+l^^' (^^^ 

where 7 is some Jordan curve about the point x = this point lies on the real axis. 

Here, we conclude the principal part of our derivation. This sequence of transformation also is need for the sequel, at least 
twice. 

Now, recall ( |34] |. Interchanging the order of summation and integration; then interchanging the order of summation (with 
proper justification) yields 



^ ED-')-- 

X e{x^\Y- 



^ / TO — /i\ /TO 



m— max(u+/i,^;) 



U \V 



where 



s = . (40) 

The next step is interchanging the outer sums order for diagonal summing as ^^^=0 X^IILo ^(""i ~ S"Lo^ X]t"=o " -''^(" + 
w, w) + X^IILi X]™=o ^('^j ^ + ''^)- Now, we have to redefine u :— u + w and w := w for the first double sum; also u :— w 
and V := V + w for the second double sum. Then, by a.^, a^, a„, and we denote a, a — wq, and /3 for the first and 

for the second double sums as shown in Table HV] Similarly, denote by fu,w{s) and by fv,w{s) the inner sums over to and 
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rewrite the result as 



a; + 1 ^ 



X 

2TTi 



where 



and 



OO oc> 
i;— 1 If— 

fuAs)= E f"r^)f™V.", ^>0, (48) 
' ^ \ II -\- 111 / \ 111 I 



7n—u-{-w-\-fi 



V w I \v + w I 



m — V'\-W 

It is easily shown that ( l48T t and ( [49] t are, in principal, similar to ( l36b . Arguing as above (see < l36b . . . < [39b ) we obtain 



and 



where F is a closed contour surrounding the point C, = z. Here we suppose that 

z — 1 

(assume z e C till the end of this section). 



(47) 



(49) 



^^i — ^^^^ — '''' 



1 /• (c-ir (f^)^c'-lr 
^^i ^^^^ 



z + l' 

Substituting ( fSOl l. and ( fSTT i for fu.w{s) and fv.w{s) in ( |47] l; then interchanging the order of summation over and integration 
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with respect to we get 



C^t) 



x + 1 



X {x + lf 



Z+1 



iC-ly{^-x)f^''-^'+^ 2"+i 



1 



-^(-irt""(:. + i) 
)j=i 

-|c<„|+/3„ 



z- 1 



27ri 



c-i 



■iu=0 



(52) 



where 



Denote by and the inner integrals as follows: 



z - 1 



and 



27ri 



ic-ir 



C + 1 



u)=0 



(c-ir 



c-i 



Note that for 



E^ 



ui=0 



dC 



dC 



< 1, 



we have convergent series in these integrals. 
Consider (|54] |. We can write 



1 



c + i 



(C-1) 



VC + 1 



Consider the integrand 



27r«/r C-^-/«(C^-l) 

Xn(C) 



1 - h 



■dC 



VuiO = 



MO 



ufc_i 
c+i 



C-^-MC^-i) 



(53) 



(54) 



(55) 
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We have denominator of the integrand V'«(C) = (C ^ ^ ^ ^(C^ ~ 1)) (C + 1) which has three zeros: 

Ci. = ^ (56) 

and Ca = —1, where lim Ci — z and G becomes infinite as h ^ 0. 

h^O 

Hence, for sufficiently small h, one may suppose that there is just one zero: (^i lying inside F, and the integrand has a 
simple pole inside the contour of integration with residue 



XuiC) 



C=Ci 



where <(C) = 1 - 2hC. 

Since /„ = Res^=<;i (77u(C)), then 



where 

Now, consider ( [55] l. We can rewrite 
Consider the integrand 



l-R \, fl-2h-R 



2h J \l + 2h-R 



R = ^/T^^4h7+Ih^. (58) 



Recall that w > 1, then we have denominator of the integrand ipviC) — C " z — h{C,'^ — 1). In this case we also have ( |56] | as 
zeros of the denominator with the same properties. Then ly ~ Res;;^^^^ {flviO) or 



l-R ^'"fl + 2h-RY 



I, ^ ^ — > y--'"-"' ^ ,59, 



By substituting ( |57] i and i59[ to (|52l l we have 



a; + 1 ^ 













2i? I 



1 

X ^(-l)«<°" /"^LtlV" T-|a«l+/3u fi-2h- RY 



u=0 



1 f^-xY z + 1 f z-ll + 2h- R 



'x + lY^ u^,^r..fl-2h-RY" 



^-IVC-ly 2i? \z + ll-2h-R, 

.|a„|+/3„ 

v=l 

Recall dTTJ, (|43] i, ( |45] l, (l42T i. (l44l i. and (|46l ) (see Table HVT l, then according to a„ and properties, we need to break these 



x^(-irt-(^) r'-'+^-l^ ^ " ) He (60) 
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series on two parts each; one finite and one infinite series. Define these series as follows: 



q 



u=0 



g-pl+P2 

9 



max 



E 

CSO 

E 



(•r, 



(•r- 



t;— max — 



g'-pi+P2 
9 



By a denote the shifted charge; then by n denote the quotient of a and q 



(T = a - pi + p2, 
a 

n 



Consider ( l60l l as four sums of integrals along the contour 7. Each of these integrals includes one of the series above. Let 
Ji, J2, J3, and J4 denote these integrals. 



J2 



where by ^i, ^2, '^3, and ^4 we denote the integrands of Ji, J2, J3, and J4 



*i(0 



£,~x\^ z + 1 fl-2h-R 



£.-x\^-lJ 2R \l + 2h~R 



[ 71—1 



E^i"' 



ti=0 

z + 1 fl-2h-R'^^ 



^-3(0 



(C-a;)(e + l)^ 2i? \l + 2h-R 

00 

X E ^2", 

inax(n.O) 

r-'^ z + 1 /z - 1 1 + 2/i - i?\ ^ 



+ 2i? \z + ll-2h~ R 

xE^i^ 



i;=l 



*4(e) = 



* /^f z + 1 /z- 1 1 + 2/i- 



X E 

D— max( — n+1,1) 



rpV 

^2 1 



where 



then (|60] | becomes 



4/i 



A-a/2l +P1 

opj.[cr/2j +P2 
+ . . (J2 + J3). 



1-i 



Note that for 



\T2\ < 1, 



(61) 



(62) 

(63) 
(64) 



(65) 
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we have the convergent geometric series in ( |6T1 ) and ( |62] i. Therefore we can rewrite these integrands in closed form. Then 
we substitute (|53] l for h; then, in turn, for z, and ( |40] i for s. Now, for the sake of simplicity and for eliminating, where 
possible, the imaginary parts of the denominators, we perform straightforward calculations and obtain 



*i(0 



2R 



Q 

X f 1 - T, 



max(n,0) 



2R 



ba{^)YH2 + il-tiit/T)'')R 



x{n,0) 



2R 



X 1 - 



Q 

max(— n,0) 



2R 



,Hl-{l-^lT'i)R 



X T. 



max( — n.O) 



where 



hi, is then defined like (l38il as follows: 



fiit/r)" P2 + P4; 



bk = 



fe+2 



2(1 -i)- 



Now ([63T l and ( |64] | become 



ri = -(t/r)^ 



i? 



2fcfe (C^-l)' 



T, = 



Pq, . . . , P4 are polynomials of degree 2 in the complex variable 

Po=^-x-{bd+ bk) {e - 1) , 
Pi = e - X - - (e - 1) , 

P2 = e - .T - - 6d) (e' - 1) , 

P3 - e - - (1 - 2r(i-^)«) + bk) {e - 1) 



b, + bk (l-2r(^-i)''))(e^ 



P4 = e - X 

Denominator functions Q, and P are defined as 

= e - ^ - (fed (1 ~ T«) + &/C (1 - ^"'')) - 1) 



R = ±jp§ - 4bdbk ie - 1) 



(66) 



(67) 



(68) 



(69) 



(70) 
(71) 



(72) 
(73) 



The functions ( |72l ) and ( |73] l define singularities of integrands ( |66] l. ( |67l ). ( |68] l. and ( |69] l. Indeed, P have 4 roots; these all roots 
are real: 



6.2 



It V1-4(V5:?+V6^) x + 4(V&;i+v^) 

2{Vb~d + Vh)^ 
lT\/l-4(Vb^-^/5^)'x + 4(^/b^- 76^)' 

2(\/E^-\/5^)' 
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Fig. 1. Contour 7 which is the circle (or any Jordan curve) about x 




Fig. 2. Dumbbell-shaped contour about the branch cut on [^3,^1] and a circle contour about ^5. 



Note, that ( |73] | becomes imaginary on [^3,^1] and [^2,^4]; therefore ( |66] |. ( l69l l have branch cuts on these segments. 
Obviously, that Q — if £^ — x. Therefore, there might be a singularity at point £,5 = x for integrands ( l66l) . . . . , ( l69l l. It may 
be simply proven that ^1, ^3, and ^5 lie inside the contour 7 as shown on Fig. [T] Thick line inside the contour denotes the 
branch cut on [^3,^1]- 

Instead of the contour 7, we can consider two contours; the first dumbbell-shaped contour whose "balls" contain, respectively, 
^3 and ^1 and the second contour •j^ as shown on Fig. |2l In this case we can write Ji, say, by the following 



'71. v"yj- 

+ / ^-1^+ / • (74) 

J-lR J III Jlx / 

We can show that the integrals along the two circles 7^ and 7^ vanish in the limit as the radius e of these circles tends to 0. 
Indeed, we can rewrite (l73T l as 

R = {bd- 6fc)V(e-ei)(e-6)(e-6)(e-^4). 

Therefore, for arbitrarily small e, we have 

R 



R 



l?-?i|<e 



l?-?3|<f 



{bd - 6fc)v/(e-6)(e-6)(e-e4) 



±{bd - 6fc)V(c-ei)(e-e2)(e-e4) 



(75) 



(76) 



Now, we estimate integrals <f^^ ^i{^)d(, and '^liO'^^ as follows. First, substituting (fTST i and (l76T l for i? in ( |66] |. we see 
that 

(&<i(^+i)r 



2(6d - 6fc)V(e-6)(C-6)(j-e4) 



are bounded above at the point ^3 and it's nearest neighborhood as well as 



2ibd - 6fc)x/(e-a)(C-6)(C-e4) 



(Me+i)r 



Q 
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are bounded above at the point and it's nearest neighborhood. Then we show that 

J'yji ^JR VS si 

/ *im<M/ -=^. 

Using parameterization of the contours 77? by ^(i?) = + ee*'', = iee^^M with e [— vr, tt] and 7^ by ^(i?) = f 3 + ee"', 
di = iee'^dd with G [0, 27r], we have 

f-ia VS. - ^1 ^ 

and we can see that <f^^ ^'i(0'^C ^"d ^i(C)'^'? vanish as e tends to zero. 

Hence, there remain two integrals <f^^ ^i{£^)d^ and ^^^^ along the line segments just above and just below the cut, 

and one integral <f^ (^66 Fig. |2]i. 

Along the upper and along the lower segments, we can rewrite 



Ji3 Jii 



This directly follows from the fact that R is the only function which becomes imaginary into these limits of integration. 
Obviously, that ( |74] i becomes 

Ji = i / ' 5 + Res(4'i(C)). (77) 

In the same way, we get J2, . . . , J4 in the form identical with ( ITTI i. 
It is not hard to show that _ 

Pq-R 

2bkie - 1) ^ 

between the limits of integration [^3,^1], where we define ip as 



r'^/^e"" (78) 



y> = -»ln( 1 (79) 



Equation ( |79l ) directly follows from the fact that R becomes imaginary on segment [Csj^i]- 

Now, recall dTOl l and dTTl i: then consider power functions T" and T^. Taken into account dTSl l. we can rewrite these functions 

as 

Tf = — {coship + ismhip), 



T2" = (-1)" j (cosn^ + isinn^), 

where n denotes either max {h, 0) or max (— n, 0). 

Since we need only imaginary parts of ( |66] |. ( |67] |. ( |68] ). and ( |69] l, then we can rewrite the sums Ji + J4 and J2 + J3 in 
as follows; 



2n Q 



cos {flip) — (1 — fir'^) sin (hp) d^ 
R\ I 



Res(*i(C)) -Res(^'4(C)) (80) 
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and 

(-1)" r-~- fbd{^~l)Y 



J2 + J3 



27r Q \ 

cos {fnp) — (1 — /i {t/rY) sin (ni^s) 



Res(*2(e))-Res(*3(C))- (81) 



To conclude the derivation of our generating function, it remains to note that 

Res(*i(O) = 0, 

= 2(1-4 ' 

Res(*s(0) = 0, 

Res(vI/4(0) - 0. 

The proof of ( l82b is straightforward but tedious and omitted here. 

Degree of radicand of R and properties of its roots define (I8OI 1 and (ISTT l as elliptic. Therefore, generating function dSTl i does 
not have a closed form. 

Thus, we just proved the following theorem. 

Theorem 1. A recursion relation (|6]l does not have a closed form. 

Note that from the residue theorem it follows that the limits ^3 and ^1 may be replaced by ^2 and ^4 respectively with 
changing signs of dSOl l and (ISTT l. 

— n n+2 



D. Generating Function of Sequence C„ C„ . . . ^ Cjj 



Consider a finite power series 

n n 

C„(y) = 5] or C„(2/) = ^ C^y^"-'^/', 



(7— — n a— — n 

a even „ odd 



where y e M. 

Here we can also consider ( l33T l instead of C,'J^. The calculations that let us got (|34] | yield now 



2ri+p 



/ m—p m—p 



m=P \ ti=0 ■u=0 
m even \ 

or 
m odd 



m-p _ \ / m-p 
2 ^* W 2 



.9(2/) , (83) 



where by (/(y) we denote the inner sum as follows: 



2n+p cT-p m-p , 



(7—m 
<7 even 



2 



cr odd 

X ( > 2 (84) 

2 

Certainly, we may use the line of reasoning that has led us from ( |27] | to ( |29] l. However, we choose another way. It is easy to 
observe that n and a form one term ^^^^ or ^^^^ (see (l33Tl). So we can repeat the steps between (|33] | and (|39] | and obtain 

n+p m—p 

5(2/) = - l)"+™-''+i(a; + 1)1^^ 

X i ii^!^^^ d£, (85) 

ZTTZ _ + ^^^^^ 
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where we redefine a and a as 

a = n — (m — p){d + 1) — pi — P2 — (u + v)q, 
a = pi + uq. 

Note that in this case we do not need the hypergeometric transformation that we have performed above. Indeed, we akeady 

2 

2n+p n 



have finite series by definition (|84| |. If we replace ^-^^^^ by to, we get 



m=P m=0 
m even 



m odd 



After interchanging the order of summation and integration in ( I83I I. we see that each of two inner sums can be identified with 
the binomial theorem. 

n-f p 



27ri (^2 - 1)(^ - 2;)n-pi-P2 + l 

X 1 - /i+ > ^-^ 

x(l-('^VV" ' il-i"—^] ) He (86) 



X + 1 J J \ \x — \ 
Now, only a finite geometric series remains in ( [86] l. Redefine ft, as a term of this series 



ft = 



(a;2-l)'^(e2_i)"-pi-P2+i 



x + ly / V Va; — 1 



Substituting '-^"j^ for ErUi ^™ •ESI' we get 



n+p 



27ri 



/7 {e-m-xr-p--p-+^'^^ 



where by ^(^) we denote the integrand as follows: 



« + " 



where 
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such that 1 — h — — ^^t+i(-^2 — jy- Now if we recall that the contour 7 was initially considered as a small circle surrounding 



the point x, then, using the residue theorem, we get 

1 
2^ 



7 

2(fe+2) lim h 



X 



5-/ 2(fe+2) 



m— 1, 



\ I ~? Fl^' otherwise, 



(88) 



where 



^ ' 1, else, 



and ^i, ^2, • • •, C2(fc+2) are roots of V'(C)- 
It can easily be checked that 



^ (f - 0. (89) 



We might prove it using the contour integration methods. But more simply is to observe that = whenever a > n. 
It remains to find the first integral in ( |87] i. Obviously, 

i^/ ^±ir 



Note that these integrals become equal to whenever n — pi — p2 + 1 < 0. 

Substituting ( l90l ). dSSl ). and ( |89] | for the first, for the second, and for the third integrals in dSTl i. we obtain the final result. 

V. Algorithms for encoding and decoding constant- weight and constant-charge binary run-length 

LIMITED sequences 

Recall that ,^ denotes a set of constant- weight or constant-charge binary run-length constrained sequences x = (xi, X2, . . . , 
Xn) of length n. Let the set ,y be ordered lexicographically. From [3] it follows that the lexicographic index of a; € is 
given by 

n 

N{x)^Y.^jW{p), (91) 

i=i 

where W{p) denotes the number of sequences in 5^ with given prefix p = (xi, X2, . . . , 0). 

The decoding algorithm, for given sequence x, find its lexicographic index A^, < iV < |c5^|. This is done by successive 
approximation method [3], [15] using ( |9T] l and W{p) as the number of sequences in .5^ with given prefix p — (xi, X2, ■ ■ ■ , 
Xj^i,0). 

By aj{p) denote the number of trailing zeros of the prefix p. By i'j(p) = i^j-i = J2i=i and by aj{p) = J2i=ii~^Y' 
denote the weight and the charge of this prefix. Since p is the prefix of x, then subsequence x = {xj,Xj+i, . . . ,a;„) is the 
rest of X, and Ij is the leading run of zeros in this subsequence. We define Ij as the complement of aj{p) in I (for leading 
run of zeros) or in k (for another run of zeros) as follows: 



/ — aj{p), = 0, 

k — aj{p), otherwise. 



Similarly, by rj we denote the trailing run of zeros in the subsequence x. If this subsequence consists of zeros, this may mean 
either the leading run of zeros (when = 0), or the trailing run of zeros. Therefore, we define rj as the complement of aj{p) 
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either in min(Z, r), or in r. In the case of nonzero x, we define rj as r, i.e., 

niin(Z, r) — aj{p), v = i^j-i = 0, 
Vj = Ir- aj{p), V = Vj_i ^ 0, 

or 

'min(Z, r) - aj{p), crn-j{p) = n - j and vj-i = 0, 
r - aj{p), an- J (p) =n- j and ^ 0, 

(^n-j{p)^n- j, 

where 

a„_,(p) = (-l)''-n'^-^.--i)-l- 
Then we can compute the number of constant-weight sequences Wi,{p) as 

K^^'j {d, k, lj,rj), Ij > and rj > 0, 



W^p) 



0, 



otherwise 



and the number of constant-charge sequences (p) as 

W,{p) = i<^n-7^''\d,k,lj,rj), Ij > and r, > 0, 
1 0, otherwise. 

Let 7 = (70, 7i, . . . , 7„) be a binary vector; components of this vector jj, < j < n indicate that x does (if — 1) or 
does not (if 7^ = 0) have the weight v = j or the charge a = 2j — n. 
Hence, the number W{p) of sequences with given prefix p be 

^^7!^W^(p), constant- weight case. 



W{p) = < 



7{a+n)/2W„{p), constant-charge case. 



a=—n 
a even 



V a odd 



We introduce the next variables: a, w, c, which correspond to aj{p), Vj-i, crj-i. 

N:=0- a:=l; w; := 0; c := 0; 
for j := 1 to n do 
Get W(p); 
if Xj = 1 then 
N — N + Wip); 
a := 1; w := w + 1; 
else 

a := a + 1; 
end if . . . else 
c:=c+(-l)'"; 
end for. 

The encoding (inverse) algorithm, for given lexicographic index N, < N < \^\, find the corresponding x. 

a := 1; w ;= 0; c := 0; 
for j := 1 to n do 
Get W{p); 
if iV > W{p) then 
N ■.= N -W{p); 
Xj—l; a := 1; w:=w + l; 
else 

Xj := 0; a := a + 1; 
end if . . . else 

c:=c+(-l)"'; 
end for. 
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TABLE V 
An Example of Single Peak Shift 



N 


X\ X2 X^ X4 Xs X6 X7 XS 


U (T 







1 1 

1 -1 -1 -1 -1-11 1 


2 

-2 


right 
shift 


1 


1 1 
1-1-1-1-11 1 1 


2 




true 
position 


2 


1 1 

1-1-1-11 1 1 1 


2 

2 


left 
shift 



VI. Further remarks 



A. 

As it was shown in Section [III the recursion relation (|6]l has an impHcit mutual nature. We can cite an example of practically 
identical coding scheme [16] in which we obtained mutual recursions in explicit form. Alternating runs in [16] were presented 
by series of zeros or ones with independent constraints. It seems now; there was a heavy construction without significant 
preference of present constant-charge code. 

B. 

Channels may cause peak shifts in dfc-constrained sequences [17]. It is considered as a more frequent error. We provide 
an example (see Table |V]) which shows that the charge distribution seems rather suitable for the peak shift analysis then the 
weight distribution. On the other hand, the weight distribution remains useful for erasure and insertion control. In more detail, 
consider the integer (or composition) representation of RLL sequences [18]. This representation is used in error detection and 
correction technique including the peak shifts correction [19], [20]. In such case an RLL sequence is parsed uniquely into 
a concatenation of phrases, each phrase beginning with one. By (f)j denote the length of jth phrase. We can state a simple 
relation between <j)j and charge aj as follows: 

4'j =J- 



where i and j are positions of consecutive ones such that i < j. Moreover, from Section and Section |V] it follows that 
Cover's enumerative scheme implies counting of (f>j. 

C. 

Since run-length constraints bound the weight distribution, we have v S [t'min, t'max], where < Vmin < i^max < n. In [6], 
Ytrehus obtained 

'0, n < min(Z, r), 

1, min(Z, r) < n < I + r + 1, 

n — I — r ^ I 



k + 1 



1, l + r + l<n 



and 



d+1 



for dklr sequences. Similarly, for the charge distribution of dk sequences we can write 

'n-2m{d+l), < n - m{k + d + 2) 
<k+l, 

|0-|max = { 2(m+ 1) 

X (A: + 1) - n, k + 1 <n- m{k + d + 2) 
<k + d + 2, 

where 



k + d + 2 
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D. 

The charge of the prefix aj (see Section HI] and Section |V|l is called the running digital sum (RDS), see [21]. The RDS has 
a finite range of values. For the RLL sequences with zero accumulated charge (dn — 0), we can see that absolute value of 
RDS does not exceed [|cr|niax/2j + 1. 

We can impose restrictions on the range of RDS values. Denote by z9i and d-z the lower and the upper bounds of this 
constrained range. The range of RDS bounded by di and i?2 is said to be the digital sum variation (DSV), see [22], [23]. 
After the authors [23], DSV constrained RLL sequences are called DCRLL sequences. 

By C^{d, k, r, -di, denote the number of DCRLL sequences beginning with one. Also by C'^{d, k, I, r, 'i?2) we denote 
the number of DCRLL sequences beginning with a leading run of zeros. Below under C^(z9i, 1^2) we consider C!^{d, k, r, -di, ^2) 
and under we similarly consider C^{d, k,l,r,'di,d2)- 

As shown in [24], calculation of the number of DCRLL sequences can be performed as concatenation of two subcodes; the 
volume of each subcode is found recurrently in explicit mutual form. 

The recurrent method for calculating the number of constant-charge RLL sequences, which we suggest in Section HIl allows 
us simply turn to calculating the number of DCRLL sequences. Indeed, we know the charge value of the prefix p at the each 
level of recursion (|6]l. Therefore we can control DSV by using the additional condition i9i,„ < (7,„ < i?2m- For this reason, 
we do not account C^™ for which this condition does not satisfy. Here under dim and '02m we consider i9i and i?2 after 
justification at the each level of recursion (|6]l. Moreover, each level of this recursion alternates the direction of charge changing; 
then it is sufficient to check either i?i,„ < (T,„ or (t,„ < d2m condition depending on the direction of charge changing. By 
substituting -di for 1^2 and vice versa, when calling C^{-di, ^2), we keep the only condition di„-i < and obtain the implicit 
mutual recursion similar to 

From initial conditions (|5]l it follows that a unique sequence of zero length and zero charge exists. In turn, from this statement 
it follows that we need an additional condition; this condition allows us to take into account the existence of Cg i?2)- Indeed, 
the sequences, which beginning with zero, have initial charge equal to 1, the sequences, which beginning with one, have initial 
charge equal to —1, and the sequence of zero length, have initial charge equal to 0. We now write this triple condition using 
the Iverson bracket notation ^ 

1, the condition is true, 
0, the condition is false. 



[condition] 



Thus we can rewrite Proposition |2] as 

and can 



Proposition 3. The numbers C'^{i)i, ^2) ond can be obtained as. 
//-[n 7^ 0] > i?2, then 



If (J = —n and the sequences begin with one, then 
If (7 ^ ~n and the sequences begin with one, then 



1, n < min(r + 1, — I?!) 
0, otherwise. 



min(n,A:+l, — i9i ) 

E C^.:-7'M2-J- d+l<n, 
j=d+i 

0, otherwise. 



If i?i > [n 7^ 0], then 

If a = n and a leading series is running, then 

Act _ 

If a ^ n and a leading series is running, then 



Q=0. 

1, n < min(/, r, 1^2), 
0, otherwise. 



min(nj,^2) 
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Evidently, the algorithms from Section |V] are suitable for encoding and decoding DCRLL sequences if Wa {p) will be 
calculated as 

'C"^-7^''\d,k,lj,r,,di,,^2j), Ij > and r, > 0, 

and 

^0, otherwise. 



where 



z9 



2j 



/ J if i^j-i is even, \ 

I [—'^2; otherwise; / 



E. 

One can find in literature some examples of application of generating functions for RLL coding. From 48 years range of 
publications [25] - [26], we cite two examples. Kolesnik and Krachkovsky [12], when deriving the Gilbert- Varshamov bound, 
estimated the volume Vr{x) of sphere in ^{n) centered on x as 

Vr(x) < mm . , 

^ ' - o<y<i yM„(0)' 

where An{y) may be taken from (l30l l. Note, that practical applications of the generating functions for enumerating RLL 
sequences often require a rational expression. Ferreira and Lin [19] and some other authors [4], [7] obtained a generating 
function for enumerating dk and dkr sequences 

t(l - f+M 
.y(t) = - I 

Assume y = 1. Then from (l28T l we immediately get the same. 



F. 

Recall (|57] |. (l58T l. and ( |59] |. One can identify /„ and /„ as generating functions of the Jacobi polynomials [14]. Then cos {fup) 
and sin [flip) in dSOl l and in dsTl l can be identified with Chebyshev polynomials in the variable ■ The technique 

which uses the theory of orthogonal polynomials is known in charge constrained coding [27]. 



G. 

We obtain the generating function ( ISTT l in the form of elliptic integrals (ISOl l and dSTT l. Further exploration may require a 
canonical form of these integrals. Using factorization, it is not so hard to reduce the first term of the integrands in dSOb and (ISTI ) 
to rational functions. In such case, we shall use polynomial expansion of cos {fnp). The next transformations, which bring 
these integrals to canonical form, may be taken from [14]. 



H. 

If we need a two-variable generating function for the case of constant-charge sequences 

c{t,y)= J2 E O'-'^^/'y^^-"^/', 

(T — Ci 71 — p 

a even n even 

or or 
a odd n odd 

then we can write 

C2 



a—c\ 
a even 



a odd 

In such case we may perform summation of ( [66l t, (|67] t, (|68] t, and < [69l t over the range [ci, C2] of even or odd numbers. In resuh 
we do not obtain any additional poles with nonzero residues. So, there remain the contours of integration depicted on Fig [2] 
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VII. Conclusion 

We have presented constant-weight and constant-charge run-length constrained binary sequences. On the base of Cover's 
enumerative technique, we have obtained recursion relations for calculating the numbers of these sequences. For investigation 
of the asymptotic behavior of these values, we have derived generating functions for enumerating such sequences. We have 
proved that generating function for enumerating constant-charge sequences does not be expressed in a closed form. So, we 
have presented the long chain of derivation steps which led us to expression for this generating function in the form of elliptic 
integrals. Also we have described two algorithms for enumerative encoding and decoding constant-weight and constant-charge 
sequences. Then, we have provided some examples of application for our results; in particular, we have extended our results 
on RDS constrained sequences. 
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